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Abstract 

> 

C"- ■ The density matrix of a spin s is fixed uniquely if the probabilites to ob- 

tain the value s upon measuring n ■ s are known for 4s(s + 1) appropriately 
^r) ' chosen directions n in space. This non-redundant set of numbers are just ex- 

pectation values of the density operator in coherent spin states, and they are 
easily obtained in an experiment carried out with a Stern-Gerlach apparatus. 
Furthermore, the experimental data can be inverted explicitly which is neces- 
sary for a parametrization of the statistical operator by the 4s(s + 1) positive 
parameters. 

c 

If an infinite ensemble of identically prepared quantum systems is available, the rules 
of quantum mechanics allow one to extract complete information about the quantum 
state. The required inversion, expressing the statistical operator entirely in terms 
of measurable quantities, is not straightforward from a mathematical point of view, 
and the data acquisition in the laboratory is challenging since it is neccessary to deal 
reliably with individual quantum systems. 

The experimental verification of a state reconstruction scheme requires high stan- 
dards in the preparation of individual quantum systems. By now, these standards 
have been met for quantum systems such as an electromagnetic wave [0], vibrating 
molecules , ions caught in a trap |J , and atoms moving freely in space after scatter- 
ing from a double slit 0. State reconstruction is reviewed in || and ||, for example, 
both from the theoretical and experimental point of view. 

State reconstruction for a quantum system with a finite number of states appears 
to be slightly easier than for particle systems living in a Hilbert space with countably 
infinite dimension. Various answers to the problem have been obtained for both mixed 
and pure states of spins with length s. Pure states in two- or three-dimensional Hilbert 
spaces (s = 1/2, 1) have been treated in a straightforward way while the general case 
is more complicated f?|, ||. Using Feynman niters, a phase sensitive version of a Stern- 
Gerlach apparatus , one can determine directly moduli and (relative) phases of the 
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individual matrix elements of the density operator [10] describing a mixed spin state. 
As shown in [O, the expectations of 4s(s + 1) linearly independent spin multipoles 
fix a unique density operator; however, no method has been indicated how to de- 
termine experimentally these values. An experimentally more realistic approach has 
been proposed in [|12|]: the density matrix is determined if 2s(4s + 1) real numbers 
are measured by using a Stern-Gerlach apparatus oriented along various directions in 
space. In ||13|| , a reduction to the minimum number of 4s (s + 1) measured probabili- 
ties has been proposed which, as shown in the Appendix, unfortunately is erroneous. 
Alternatively, a tomographic approach adapted to finite-dimensional Hilbert spaces 



allows for experimental reconstruction of quantum states [14 . As is common for meth- 



ods involving Wigner functions, the information to be extracted from experiments is 



redundant. For an application of this approach, see |L5|], where the determination of 
a single quantized cavity mode is treated. 

In this letter, a simple and constructive reconstruction scheme for the mixed state 
of a spin s is presented in accordance with the following two natural constraints taken 
from 



(1) the measurements are performed with a standard Stern-Gerlach apparatus only; 

(2) no redundant information is acquired. 

These conditions guarantee that a standard experimental setting is sufficient to per- 
form the reconstruction of mixed spin states in the most economical way. 

The states of a spin of magnitude s belong to a Hilbert space Ti s of complex di- 
mension (2s + l), carrying an irreducible representation of the group SU(2). The com- 
ponents of the spin operator S = hs with standard commutation relations [s x , s y ] = 
is z , . . ., generate rotations about the corresponding axes. The standard basis of the 
space 7i s is given by the eigenvectors of the z component s z = n z ■ s of the spin, which 
are denoted by n z ), — s < p < s. The phases of the states are fixed by the transfor- 
mation under the anti-unitary time reversal operator T: T\p, n z ) = (— 1) S ~ M | — /i, n z ), 
and the ladder operators s± = s x ±is y act as usual in this basis: 



s±\fi,n z ) = yjs(s + 1) - ± ± l,n g ) . (1) 

The complexified algebra A s of observables in the space Ti, s has complex dimension 
(2s + l) 2 . It consists of all polynomials in the operators s^, s y and s z with complex 
coefficients and of degree 2s at most. A monomial of a degree higher than 2s can 
always be expressed as a linear combination of monomials of lower degree. 
Consider the eigenstates of the operator n ■ s, 

n ■ s [//, ft) — /t |/t, ft) , —s<fi<s, (2) 

where the unit vector n = (sin 9 cos ip, sin 9 sin ip, cos 8), < 9 < 7r,0 < tp < 2tt, 
defines a direction in space. Given a state with density matrix p, the probability 
P)i{n) to measure the value \i with a Stern-Gerlach apparatus oriented along n is 
given by 

p^(n) = Tr [pPf,(n) } = (fi, n\ p n) , (3) 
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with the projector P^(n) = \p, n)(p, n\. For the state with maximal weight, p = s, the 
probability p s (n) coincides with the expectation value of p in a (standard)^] coherent 
state [fP7|1 , 

|s, n) = exp[— i 9 rh((p) ■ s] |s, n 2 ) , (4) 

where m(</>) = (— sin ip, cos 0). In other words, the state |n) is the result of rotating 
the state \s, n z ) about the axis m((p) in the xy plane by an angle 9. It is convenient to 
combine (9, ip) into a single complex variable, z = tan(6 l /2) exp[i<£>]. This corresponds 
to a stereographic projection of the surface of the sphere to the complex plane. In 
terms of z, a coherent state has the expansion [18] 



According to (|3|) the mean value of p in this state is 

p " (n)= „,£A-"'i I s -») ^wr^" (6) 

where z denotes the complex conjugate of z. Note that Eq. (|j) defines a linear 
relation between the quantities p s (n) and the unknowns p^ ^. Due to the hermiticity 
of the density matrix, p^ ^ = p„„>, it contains (2s + l) 2 free real parameters (the 
normalization condition Tr [p] — 1 will be suppressed for the moment). Therefore, 
the probabilities p s (n) must be known for at least (2s+l) 2 points z\, A = 1, 2, . . . , (2s+ 
l) 2 . Each of the points z\ in the complex plane corresponds to a point on the sphere 
or, equivalently, to one spatial direction n A = n{9\, <p\). 

If the points z\ are chosen appropriately, the (2s + l) 2 measurable numbers p s {n\) 
contain all the information needed to determine the matrix elements p^ ^ and thus 
the quantum state. To show this, let us rescale both the measured probabilities, 

Px = (l + \z x \ 2 ) 2s p s (n x ), (7) 

and the matrix elements of the statistical operator, 

. _ ( 2s \" ( 2s X' 2 , R , 

Pk'k — \ fc> J \ k j Ps-k's-k' \°) 

which, for convenience, have been relabelled with nonnegative integers k = s — p and 
k' = s — p', k, k' = 0, 1, . . . , 2s. Writing down Eq. (^) for (2s + l) 2 different (as yet 
unspecified) directions fix, one obtains 

p x = J2 ^4 Pk'k, A = l,2,...,(2s + 1) 2 . (9) 

k,k'=0 



1 Other families of coherent states are obtained if a 'fiducial' state different from \s,n z ) is used 
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It it not obvious how to invert directly the (2s + l) 2 x (2s + l) 2 matrix Nu^/n = ~z k \z\, 
which would provide an immediate solution of the problem. By a Fourier transform, 
however, one can divide the (2s + l) 2 coupled equations into smaller sets of equations 
each of which contains (2s + 1) unknowns. For integer spin, (s + 1) such sets will 
emerge while for half integer spin their number is (s + 3/2). Select (2s + l) 2 directions 
n\ = n qr with 

z\ = V = R q expfivv] > < q,r <2s , (10) 
with positive numbers R q > 0, R q ^ R q i if q ^ q', and 

2n 1 

<p qr = (r + qA), 0<A< . (11) 

yqr 2s + V H ! ~ ~ 2s + 1 V ; 

The numbers R q define (2s + 1) circles in the complex plane which correspond to 
(2s + 1) circles on the sphere about the z axis. The values of the angles ip qr define 
(2s + 1) directions equidistant on each circle. An angle A ^ provides a shift of the 
directions on one circle relative to those on the others. It turns out that a nonzero A 
is only necessary if one deals with the fermionic problem of state reconstruction (half 
integer spin) while it can be dropped in the bosonic case (integer spin). 
Using 

1 2s+l 

— — - exp[i(m + k-k')ip qr } = 8 k > k+m + exp[i27rgA]<5 fc / fc+m _ (2s+ i) , < m < 2s , 
Zs + L r=1 

(12) 



Eq. (|) turns into 

P,m = E K +m Pk+mk + exp[*27rgA] £ R 2k+m- { 2s+i) h+m _ [2s+l) fc ; (13) 

fc=0 fc=2s-m+l 

where the shorthand 

I 2s+l 

Vqm = exp[im(f qr }p qr (14) 

is + 1 r=1 

has been introduced. A matrix notation will be useful here. Collect the unknowns 
associated with a fixed value of m into a vector with (2s + 1) components, 

Pm — {PmOj Pl+ml; • • • > P2s2s-m] P02s-m+l, ■ ■ ■ , Pm-12s) , < 111 < 2s , (15) 

(with no entry on the right of the semicolon if m = 0) and similarly the data p qm into 
vectors 

Pm = (plm,P2m, ■ ■ ■ ,P2s+lm) , < 1TL < 2s . (16) 

Then, the relations ( |13|) take the form 

p m = M m p m , 0<m<2s, (17) 

where the elements (M m ) q k of the (2s+l) matrices M m with dimension (2s+l) x (2s+l) 
can be read off from ([13]). For each value of m, Eq. ( p[) is a closed set of equations for 
the unknowns p m . If m equals zero, there are (2s + l) real equations for the (2s + l) real 
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diagonal elements po of the statistical operator. If m equals 1, one has (2s + 1) linear 
complex equations for (2s+l) complex numbers v%. The remaining sets of equations for 
m = 2, 3, . . . , 2s have the same structure. However, not all of them are independent. 
In the bosonic case, the sets come in s pairs: taking the complex conjugate of Eq . (|17D 
with label m , one obtains the equation with index (2s — m ). Hence, the total number 
of real independent equations is at most 1 ■ (2s + 1) + s ■ 2(2s + 1) = (2s + l) 2 which 
coincides exactly with the number of real free parameters of the density matrix p. In 
the fermionic case, there is again one set (m = 0) of real equations for the (2s + 1) 
diagonal elements of the statistical operator, plus (s — 1/2) paired sets of complex 
equations, plus one 'self-conjugate' set of equations with only (2s + 1) real unknowns. 
This gives a total of 1 ■ (2s + 1) + (s - 1/2) ■ 2(2s + 1) + 1 ■ (2s + 1) = (2s + l) 2 , as 
before. Thus, it remains to show that the matrices M m are invertible for the relevant 
values of m. 
The choice 



R„ 



S—q 



r > 



(18) 



turns each M m into a matrix of type Vandermonde. The bosonic case (A = 0) is 
particularly simple: defining 



r 



2k+m 



one obtains: 



2fc+m-(2s+l) 



M, 



if < k < 2s - m , 
if 2s - m + 1 < k < 2s 

■■• (r2s) s \ 



(19) 



(20) 



V (r )- s ••• (r 2s y s ) 
The determinant of a Vandermonde matrix is known explicitly []T9 , leading here to: 



det M r 




n (r* 

0<k'<k<2s 



(21) 



and it obviously vanishes if and only if two numbers and rv, say, are equal. However, 
r being nonzero, one has r k /r k = r 2<yk ~ k ) ^ if k ^ k', and the inverse matrices M" 1 
do exist for this choice of (2s + l) 2 directions n qr . An explicit form of the inverse of 
matrices such as M m has been worked out in 

The reconstruction of pure states in systems with half-integer spin s proceeds 
similarly. However, a nonzero shift is necessary: A = l/(2s + 1), say. In this case, 
one obtains again matrices M m (0 < m < s + 1/2) of type Vandermonde: for < 
k < 2s — m the entries simply read {N\ m )qk — { r k) s ~ 9 , while the remaining entries are 
given by 



[M m ) qk = (r fc exp[-i27rA]) s q exp[i2nsA] , 2s-m + l<k<2s. 



(22) 



These extra phases exp[— z27rA] distinguish lines which otherwise would be identical, 
a coincidence which does not occur for integer values of s. Therefore, all the matrices 
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M m (0 < m < s + 1/2) are invertible, too. For a spin s = 1/2, the directions n 00 and 
Uqx are located in the yz plane, while n w and nn are in the xz plane, providing thus 
four independent numbers to determine the (unnormalized) density operator. With a 
zero shift A, one would obtain information related to the yz plane only which is not 
sufficient for reconstruction. 

In summary, it is possible to reconstruct the density matrix p from (2s + l) 2 
probabilities p s (n qr ), 1 < q,r < 2s + 1 along judiciously chosen directions n qr using 
a Fourier transform and standard linear algebra. The Fourier transform reduces the 
original problem of inverting a (2s + 1) 2 x (2s + 1) 2 matrix to that of inverting a number 
of (2s + 1) x (2s + 1) matrices of Vandermonde type. The method presented here is 
an optimal tool for reconstructing a spin state in the sense that a minimal number of 
data has to be acquired using nothing but a standard Stern-Gerlach apparatus. 
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Appendix 

In ||13|| , four different approaches have been proposed to determine the coefficients 
p* m , defined by expanding the statistical operator p for a spin s according to p = 
l/(2s + l)J2im PtmKim, where the Ki m , I = 0, . . . , 2s, —I < m < I, are a set of 
(2s + l) 2 orthogonal multipole operators [pOR . The reconstruction of the statistical 
operator p is based throughout on the inversion of the relation 

, 1 ; 

iMvO = (^t) 2 E^(^)pL, (23) 

where the functions Yi m (8, (p) are standard spherical harmonics. The functions on the 
left-hand-side of this equation are linear combinations of measurable quantities,^ 



n,(0, <p) = V2s~TT £ (-i)^(s/i, s - ti\io) Pft (e, <p) , (24) 

fl= — S 

where (sp, sp'\lm) denotes a Clebsch-Gordan coefficient, and p^O, if) is the probabil- 
ity to measure the value p if the orientation of the Stern-Gerlach apparatus defines 
n(8, ip) as axis of quantization. 

Unfortunately, the fourth method to invert Eq. ( p3"D is erroneous.^ It has been 
proposed to measure the probabilities at a fixed angle 9j = 9m, and angles ipk = 
k2n/(2s + 1), k = 0,1,..., 2s, corresponding to (2s + 1) directions located on a 
cone about the z axis. Contrary to what is stated just before Eq. (22) in fl3| , no 



2 The sums in Eqs. (11) and (13) in [[l3| both go from — s to s, not from — I to I. 
3 The authors would like to thank R. F. Werner for pointing this out. 
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orthogonality relation is available for the the angles ipk] instead one has 



Y 2s 

— — — exp[i(m - m')ip k ] = 5 mm > + <5 mm '+( 2 s+i) + 5 mm '-(2 S +i) (25) 

ZS + 1 fc=0 

since —4s < m — m! < 4s (as follows from —2s < m,m' < 2s and not — s < m,m' < s, 
as stated incorrectly), and the inversion of ( p3|) becomes impossible. Therefore, know- 
ing the values Ui(9m, fk) is not sufficient to determine the coefficients p* m unambigu- 
ously. One way out of this difficulty is to measure the probabilities along (4s + 1) 
(instead of (2s + 1)) directions distributed homogeneously on a cone. Then, relation 
is replaced by 

-y 4s 

— — - £ exp[z(m - m')A;27r/(4s + 1)] = S mm > , (26) 

and the intended inversion becomes possible. However, the number of measured pa- 
rameters has been increased considerably. Effectively, one obtains a method of state 



reconstruction which is equivalent to the approach developed in [12] 
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